Fluorescence into flat and structured radiation continua: An atomic density matrix 

without a master equation. 
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We investigate an atomic A-system with one transition coupled to a laser field and a flat continuum 
of vacuum modes and the other transition coupled to field modes near the edge of a photonic band 
gap. The system requires simultaneous treatment of Markovian and non-Markovian dissipation 
processes, but the photonic band gap-continuum can not be eliminated within a density matrix 
treatment. Instead we propose a formalism based on Monte-Carlo wavefunctions, and we present 
results relevant to the experimental characterization of a structured continuum. 
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With the advent of photonic band gap (PBG) materials 
and dispersive media, the mode structure of the electro- 
magnetic field can be tailored in a controllable fashion 
providing for instance band gaps or defect modes of var- 
ious forms (l]J^. The rapidly varying mode structure in 
the radiation reservoir invalidates the Born-Markov ap- 
proximations normally employed for a simple quantum 
system like an atom when this is located inside a PBG- 
material with transition frequency near the edge of the 
gap. The reservoir degrees of freedom are thus not easily 
eliminated to derive a master equation for the reduced 
system dynamics. The main body of theoretical 

works on atomic interactions within PBG materials has 
therefore addressed the unitary dynamics in terms of the 
complete atom(s)-|- field wavefunctions. 
In this paper, we address a A-system with one laser- 
driven transition experiencing a flat vacuum without 
structure and the frequency of the other transition near 
the edge of a PBG. A unitary wavefunction dynamics is 
incompatible with the treatment of the atomic fluores- 
cence on the "free" -space transition, and we must seek 
a way to apply the simple Markovian properties of this 
process in the solution of the complete problem. A formu- 
lation in terms of Monte-Carlo wavefunctions (MCWF) 
turns out to be particularly useful for this purpose. This 
method then also suggests itself as a means of solv- 
ing other problems emerging in the overlapping domain 
of quantum optics, semiconductors and nano-structures 
where dissipation of Markovian and non-Markovian char- 
acter may co-exist. Furthermore, our work establishes an 
application of the MCWF treatment which goes beyond 
its conventional correspondence with Born-Markov mas- 
ter equations; at no point does such a master equation 
appear in this work. 

The A-system is interesting from an experimental point 
of view since atoms may be present in their ground state 
in the dielectric host, and the dynamics of the interac- 
tion with the field modes in the vicinity of the gap may 
be studied when the laser excitation on the "free" -space 
transition is turned on. In partial analogy with the shelv- 
ing scheme technique we note that the fluorescence signal 



on the "free" -space transition may serve as to probe de- 
tails of the interaction between the atom and the field 
modes in the PBG material. 

In a PBG one finds a modified dispersion relation for 
the photons in the radiation reservoir. The methods pre- 
sented in the following apply for any mode structure and 
dispersion relation, but to illustrate the method we em- 
ploy the isotropic model introduced by John and Wang 
. Their dispersion relation for a periodic array of di- 
electric scatterers of radius a and index of refaction n, 
reads (with a separation of 5 = 2an between the scatter- 
ers): 
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which leads to a gap centered at the frequency loq 
. With n = 1.082 the gap width Auj is O.OSwg and the 
upper band edge frequency LUe is given by coe = loq + 
iAcjo — 1.025cjo. In the vicinity of the band edge we 
have (effective mass approximation) 



(1) 



where A ~ —2ac/sm(AnaLUe/c) and ko = tt/ (2a(n + 1)). 

We consider a three-level atom with two lower levels \a) 
and |c) coupled by the electric dipole coupling to a com- 
mon excited level \b), see fig. § On the I a) ^ \b) tran- 
sition we apply a laser field, and the atom may decay by 
spontaneous emission due to the coupling to a flat radi- 
ation reservoir. The transition \b) — > |c) is accompanied 
by the emission of a photon with frequency in the vicinity 
of the photonic band gap edge, and this atomic transition 
is significantly modified by the presence of the dielectric 
host. 

Neglecting the zero-point energies of the field modes, 
and setting the atomic energy levels to the values 0, huJb 
and hcuc respectively, we write the Hamiltonian for the 
system (Ti = 1) 
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H = LUbO-bb + LUcO-cc + ^ (^\aia\ + ^ i^xb\bx + V (2) 



where the interaction term in the rotating wave approx- 
imation is given by 

V = i^9\ia\(^ab - ax<Jba) + i^^9\{b\cTcb - b\abc) 

X X 

+igL{<Jbae~'^'^' - aabe'^"-') (3) 

where aij denote atomic dyadic operators \i){j\ with 
i,j (z {a, 6,c}; a\,b\ are the field annihilation operators 
of the flat vacuum and PBG vacuum, respectively, and 
the laser field is represented by a semiclassical c-number 
field. We assume that the coupling to the flat continuum 
may be treated by perturbation theory in the usual way, 
i.e. an energy shift (Lamb shift) and a decay rate 7 may 
be attributed to the excited state \h). The Lamb shift 
is assumed incorporated in the atomic energy tob in eq. 

and the decay rate describes an incoherent transi- 
tion mechanism by which atoms in the excited state \b) 
decay to the ground state \a). We shall incorporate the 
decay mechanism by an effective non-hermitian Hamilto- 
nian iJgfj- First, we identify the wavefunction evolution 
governed by this Hamiltonian, and next, by appealing to 
the Monte Carlo wavefunction formalism we shall obtain 
the exact evolution of the atomic system. 

We apply the resolvent operator deflned as |^ 
G{z) ^ 



z- H. 



eff 

where z is a complex Laplace variable and is the 
effective non-hcrmitian Hamiltonian of the system given 
by eqs. (||) and (^), but with the sums over a\ operators 
suppressed, and the replacement ujb ^ tUb — *7/2. 
With the system initially in state a, the resolvent opera- 
tor equations read 



(Z - 0)Gaa{z) = 1 + VabGbaiz + UJl) 
{z ~ UJx)Gc^a{z) = Vc^bGba{z) 



(4) 



{z-UJb + -Z-)Gba{z) = VbaGaa{z - UJl) + Vbc>^Gc^a{z) 



where Vab = gL and the amplitudes Gc^aiz) pertain to 
the PBG-continuum states |c) |1a)- Using eq. (|l|) in 
the summation over continuum modes and turning the 
summation into an integral, we get 



(z - Wf, + i-f/2)Gba{z) = VbaGaaiz - UJ l) - 



iGGbajz) 
\/z - CJe ■ 



where the effective dipole coupling to the mode structure 
is given by C = d^fcoLL'e/(47reoVC4) with d the atomic 
dipole moment on the b c transition. 
Solving these coupled, algebraic equations for Gaa and 
Gba, we find 



Gaaiz) 



Gbaiz) 



{z -LUb + ihl) + iG/^/z - UJe 



{z - ujl)[{z -u}b + ihl) + iGj^z ~ u}^\ - iKbp 



Vba 



(z - UJl)[{z -UJb+ 157) + iC/ y/z - We J - \ VbaV 



The dynamics of the system is obtained by inverting the 
amplitudes to time domain by means of the inversion 
integral for the time evolution operator. 



dzG{z)e 



where e is an infinitesimal small positive quantity. 
Due to the high order of the polynomial of z in the de- 
nominator and the presence of the square root terms in 
and (^), it is not easy to apply the residue theorem 
and to obtain the amplitudes in time domain analyti- 
cally. Instead we compute the two inversion integrals 
numerically. This integration is straightforward, and the 
calculation yields for example the populations of the ini- 
tial ground state \a) and of the excited atomic state \b) as 
function of time, ^O(i) = \Uaa{t)\'' , <(t) = \Uba{t)\'' . We 
keep track of the norm P{t) of the wavefunction, noting 
that it changes only due to the imaginary part of the ex- 
cited state energy, and hence ^\ioss — — 77i'b(0; which 
in integrated form reads 



(7) 



In fig. ^ we show an example of the relevant time depen- 
dent quantities P{t), 7r°(i) and 7r^(i). From the figure, 
it is evident that the populations in the states \a) and \b) 
approach zero after a transient evolution. There is, how- 
ever, a substantial part of the population (P(oo) ~ 20%) 
which is not lost by fluorescence on the free-space tran- 
sition. This population is transferred to the atomic state 
|c) associated with the PBG-continuum, and at any time 
we have 7r,"(i) - P{t) - 7r^(t) - 7r,"(t). 

The spontaneous decay on the |6) — > \a) transition was 
treated only as a loss mechanism for the excited state am- 
plitude, but the atoms are incoherently fed back in the 
ground state |a), and from here they are re-excited by 
the laser. In a density matrix formulation, when tracing 
over the resulting different photon number states of the 
fiat reservoir, the elimination of the PBG-modes would 
be exceedingly difficult, if possible at all. 
It has been shown that dissipative problems in quan- 
tum mechanics may be solved by stochastic wavefunction 
equations as an alternative to master equations [|o|-^. 
In the "quantum jump" scheme, one propagates state 
vectors according to a non-hermitian Hamiltonian, and 
at certain instants of time, chosen according to a ran- 
dom process, this propagation is interrupted by quantum 
jump projections of the state vectors (see Q for a re- 
cent review). In the formulations of the method so far, 
the continuously propagated state vector is described as 
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the solution of a Schrodinger equation, but, the atomic 
populations 7r° (t) identified after elimination of the PBG 
reservoir above may be applied just as well for the con- 
struction of the atomic density matrix. 
The function P{t) is the norm of the no-jump state vec- 
tor jl^ and consequently the probability that a photon 
has not been registered in the flat reservoir at time t. 
The ensemble averaged populations can be found by so- 
lution of integral equations: The population of an atomic 
state Tfi(t) is a sum of a term representing the popula- 
tion given that the atom has not decayed and a term 
representing the population given that the latest jump 
occured at time t', 

(t) = {t) + f dt'jW, {t')n^ {t-t'). (8) 

Note that 7r^"(t) = n^{t)/P{t) ■ P{t) provides the given 
normalized population with the appropiate no-jump 
weight-factor. Eq.(||) must be solved for 7fb(t) first, e.g. 
by a Laplace transform: 7f{,(z) = 7r°(z)/(l — 77r^(z)), and 
one may subsequently obtain the other populations (see 
also 11). 

The populations can of course also be found by simula- 
tions. In a single trajectory one considers the normal- 
ized populations 7ri(i) = ■K^{t)/P{t) (and other density 
matrix elements if necessary), until a jump occurs when 
P{t) equals a random number e chosen uniformly on the 
interval betwen zero and unity. The quantum jump puts 
the atom in the state |a), and from here the evolution 
starts over again. In fig. 0, we plot Tfa{t) and 7ff,(i) ob- 
tained by an average of 1(P stochastic wavefunctions for 
the same parameters as used in fig. ^. 

Let us comment on the atomic dynamics obtained in 
fig. H After an initial transient evolution, the popu- 
lations T:a,b{t) approach zero in a non-exponential way. 
The fiuorescence signal on the \h) \a) transition thus 
vanishes as opposed to the case of a two-level atom in 
free space. In the simulations we note that no jump will 
occur if the random number e is smaller than P(oo). 
When this value is non-zero, each realization only ex- 
hibits a limited number of jumps since eventually the 
value chosen for e will be smaller than P{oo). The prob- 
ability of having exactly k photon emissions (jumps) 
in a given simulated trajectory is (1 — P(oo))'^P(oo), 
and the mean number of photons emitted per atom is 
P{oo)~^ — 1. The amplitude Wc;^a(0 is the inverse Laplace 
transform of Gc^a = Vc^bGba/ {z — uj\). In the long 
time-limit only the pole z = ivx contributes and thus 
Kcait) = Vc,bGbai^x)e-'^^* and since P{oo) = 7r°{(x) 
we find 

P{^) = Y,\V,,tGtaMf 
A 

The summation can be turned into an integral, which is 
calculated numerically. 



For a A-system in free space with a branching of the de- 
cay from the upper state, there will also be a finite num- 
ber of fluorescence photons emitted on the laser driven 
transition, corresponding to P(oo) = l' /{l + 7') with 
7,7' being the decay rates of state |6) to the states \a) 
and |c) respectively, and the total number of fluorescence 
photons on the \b) — > \a) transition is thus independent 
of the parameters of the driving field. This is different 
in our case, as seen in fig. |[ where P(oo) and the mean 
number of flourescence photons emitted on the free space 
transition are plotted as function of the laser detuning 
from the PBG edge for different choices of the laser cou- 
pling. For a rather weak laser coupling, the transition to 
the atomic state |c) is a Raman-process which is strongly 
suppressed when the laser is tuned below the band gap 
edge since there are then no resonant PBG-modes for 
the Stokes photon. A stronger laser coupling leads to 
an Autler-Townes splitting of level \h) and population is 
then transferred to the PBG-continuum by a higher or- 
der process, removing the step-like character of P{oo). 
The fluorescence signal may hence probe details of the 
PBG structure. 

In conclusion we have demonstrated a technique for the 
solution of a problem, for which a Born-Markov master 
equation does not exist. We have in a parallel study de- 
rived a reduced master equation for the A-system apply- 
ing only the Born-approximation. The numerical solution 
of the resulting non-Markovian master equation yields 
a very poor agreement with the exact results presented 
here thus invalidating the use of the Born-approximation. 
These results will be presented elsewhere. The specific 
form of the structured continuum is not essential for our 
approach, but it is important that only one photon states 
of the PBG-continuum appear (a possible slow decay 
from state |c) back to \a) can only be treated if we may 
assume that the photon in the PBG continuum escapes 
before the atom is reexcited to level \b)). The simula- 
tions and the analytical expression (^) are simplified by 
the fact that all jumps put the atom in the same state. 
Our formalism, however, is perfectly capable of treating 
more general systems with branching of the decay from 
state 1 5) to multiple states \aj). This implies that the 
index a is replaced by the set of indices aj with the corre- 
sponding enlargement of the set of equations (^) and (||). 
The no-jump evolution and the associated delay function, 
following a jump to a given level Uj , are then readily com- 
puted. 

The MCWF technique has been applied to non- 
Markovian problems through the solution of Markovian 
master equations for enlarged model systems ||l^. Our 
situation, however, is different, since without ever hav- 
ing a master equation we have, by recourse to conditioned 
wavefunction dynamics, been able to obtain the atomic 
density matrix. We anticipate that by combination of 
the ideas in ref. [0 and in this paper, a wide class of 
non-Markovian problems may become tractable. 
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FIG. 2. The populations ^^{t) (dashed line), 7r°(t) 
(solid line) and the norm P{t) (dotted line) are plot- 
ted as functions of time. The parameters chosen are: 




FIG. 3. The populations Tiait) (dashed line) and libit) 
(solid line) are plotted as functions of time. The parame- 
ters chosen are the same as in fig. 2. The curves are averaged 
over lO** trajectories. 
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FIG. 1. Level scheme 
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FIG. 4. Solid line is -P(oo) for Vab = 7/2 and dashed line 
the number of fluorescence photons. The dot-dashed line is 
P(oo) for Vab = 37 and the long-dashed line the number of 
fluorescence photons. The parameters chosen: C'^^^ = 7, 

IjJb = We 
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